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We report a detailed experimental study of the structure and dynamics of glassy states in hard 
spheres with short-range attraction. The system is a suspension of nearly-hard-sphere colloidal 
particles and non-adsorbing linear polymer which induces a depletion attraction between the par- 
ticles. Observation of crystallization reveals a re-entrant glass transition. Static light scattering 
shows a continuous change in the static structure factors upon increasing attraction. Dynamic light 
scattering results, which cover 11 orders of magnitude in time, are consistent with the existence 
of two distinct kinds of glasses, those dominated by inter-particle repulsion and caging, and those 
dominated by attraction. Samples close to the 'A3 point' predicted by mode coupling theory for 
such systems show very slow, logarithmic dynamics. 



I. INTRODUCTION 

Glassy states are found in many systems [1]. How- 
ever, miderstanding the glass transition is still a major 
challenge for statistical and condensed-matter physics. 
Simple and well-characterized models hold an important 
place in this field. The glass transition in the simplest 
model colloid, a suspension of hard spheres, has been 
studied in detail for more than a decade [2-5] . The phase 
behaviour of a system of N hard-spheres of radius R in 
volume V is determined by a single variable, the den- 
sity or volume fraction — {A/i)TTR^N /V . Increasing 
(f) drives the system from a stable fluid to a fluid-crystal 
coexistence, and then a fully crystallized phase [2], which 
should be the thermodynamically favorable phase up to 
(j) = 7r/3\/2 w 0.74. However, at > « 0.58, hard 
spheres fail to crystallize [2, 3]. This is usually inter- 
preted as a glass transition due to the caging of particles 
by each other. The most successful theoretical account 
given of this transition to date is from mode coupling 
theory (MCT) [3, 4]. Within this framework, the cou- 
pling between different density fluctuation modes drives 
the system into a dynamically arrested state [5] . 

More recently, the focus of attention has moved on to 
hard spheres with a short-range attraction. Besides be- 
ing a good model for understanding the fundamentals of 
the glass transition, such 'sticky hard spheres' are also 
ubiquitous in applications. Attraction in hard-spheres 
can be realized experimentally by adding non-adsorbing 
polymers to colloids. The center of mass of a polymer 
coil of radius of gyration Vg is excluded from a zone of 
width ~ Tg from the surface of each colloid. When two 
colloids come close enough to each other so that their 
polymer-excluded regions overlap, the imbalance in poly- 
mer osmotic pressure pushes them together. This effec- 
tive 'depletion' attraction is well described [6] by the 
Asakura-Oosawa form [7]. Its dimensionless range can 
be estimated by the ratio ^ ^gl^-, while its strength 
is governed by the concentration of polymer coils in the 
free volume available to them, Cp°''. The free volume 
depends on the exact structure of the suspension, and 
is not an easily obtained experimental parameter. How- 
ever, the concentration of polymer in the whole system. 



Cp, can be used as an alternative variable to describe the 
composition of the system [8] . 

The presence of a short-range attraction in hard 
spheres widens the equilibrium fluid-crystal coexistence 
region in the phase diagram [8, 9], and introduces (non- 
equilibrium) gels at low volume fractions [10, 11] and a re- 
entrant glass transition at high volume fractions [12-22]. 
In this paper we report a comprehensive study of struc- 
ture and dynamics in the vicinity of this re-entrant glass 
transition in a model colloid-polymer mixture: sterically- 
stabilized polymethylmethacrylate particles with added 
linear polystyrene [23]. 

We used light scattering to study the structure of col- 
loids by measuring the static structure factor (SSF), S{q), 
which is effectively the Fourier transform of the pair cor- 
relation function. Dynamic light scattering was used to 
obtain the normalized collective dynamic structure fac- 
tor (DSF), /(9,t), which measures the time correlation 
of particle density fluctuations at wave vector q after de- 
lay time T. Our results agree in broad outline with pre- 
vious experimental studies and the trends predicted by 
MCT, while the detailed nature of our study and the very 
wide time window of our measured DSFs (11 orders of 
magnitude) together shed new light on the nature of the 
re-entrant glass transition in sticky hard spheres. Pre- 
liminary reports of some of these data have been given 
before [17, 24]. 

II. EXPERIMENTAL METHODS 

A. Sample preparation 

The colloidal particles used in this study were poly- 
methylmethacrylate (PMMA) spheres sterically stabi- 
lized by chemically-grafted poly-(12-hydroxystearic acid) 
(PHSA) dispersed in cis-decalin [25]. The solvated 
PHSA, approximately lOnm thick, produces a nearly- 
hard-sphere interaction between the colloidal particles 
[26]. The particle radius, R — 202nm, was determined 
from the lattice spacing of the crystal phase at fluid- 
crystal coexistence, taking the crystal to be at = 0.545. 
Particle size polydispersity was 0.069, measured from the 
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apparent angle dependence of the diffusion coefficient in 

a dilute suspension [27]. 

The colloidal volume fraction was calibrated by mea- 
suring the amount of crystal phase in the coexistence 
region and taking the fluid and crystal volume fraction 
to be at 0.494 and 0.545 respectively. Samples of the 
volume-fraction-calibrated stock were also weighed and 
dried in a vacuum oven to determine the effective den- 
sity of the particles, which was then used in subsequent 
calculations of the volmiic fraction of samples prepared 
from the stock. The stock volume fraction was also re- 
calibrated from time to time by drying and weighing. 

To induce attraction between the colloids, we 
added non- adsorbing linear polystyrene. This well- 
characterized model colloid-polymer mixture has been 
studied extensively over the last decade [23]. The poly- 
mer used was purchased from Polymer Laboratories and 
had a molecular weight of = 370000 daltons. fts 
radius of gyration in cis-decalin at 20°C was calcu- 
lated from the data in [28] to be Vg = 17.8 nm. Thus 
the dimensionless range of the depletion attraction is 
^ ^ rg/R = 0.09. 

Colloid-polymer mixture samples were prepared by 

mixing stocks of colloids and polymers with known con- 
centration and solvent by weight. Sample volumes were 
about 1 cm'^. Knowing the density of each species per- 
mits calculation of the final composition. 

The main uncertainty in sample composition comes 
from a systematic uncertainty in the calibrated volume 
fraction of the colloid stock. This is because the vol- 
ume fractions of coexisting fluid and crystal phases for 
slightly polydisperse hard-spheres are slightly different 
from those in a monodisperse colloid, but the exact val- 
ues are uncertain [29, 30]. However, all samples were 
prepared from the same stocks of colloids and polymer 
solutions, or stocks calibrated against each other. Some 
samples were also derived from others in a controlled way 
(see next paragraph). Therefore despite some systematic 
uncertainties in the estimation of absolute volume frac- 
tions due to polydispersity, the uncertainties in sample 
compositions relative to each other were mostly from ran- 
dom errors in weighing. These uncertainties are below a 
percent in the worst case and are insignificant in this 
work. 

Samples were tumbled for prolonged periods of time to 
ensure proper mixing of the components. After homog- 
enizing, a small amount of each sample was transferred 
to 3 mm-inner diameter glass tubes and sealed for light 
scattering experiments. The rest of the sample was then 
left undisturbed for visual observation of any phase tran- 
sitions until sedimentation appeared. Then some sam- 
ples may be diluted with solvent while others were left 
opened for solvent to evaporate before re-homogenizing, 
fn this way, a sequence of samples, some very close in 
composition, could be prepared. 



B. Light scattering methods 

The slight difference in the refractive indices of PMMA 

and cis-decalin (1.49 and 1.48 respectively) was enough 
to render all of our samples turbid (transmission coeffi- 
cients « 20-40%). We therefore used two-color light scat- 
tering to extract the singly-scattered component. The 
detailed experimental arrangement and data analysis for 
this method have been described elsewhere [31]. Here we 
just summarize relevant procedures. 

Two lasers of different wavelengths, blue (Ag = 488 
nm) and green (Ag = 514.5 nm), and two detectors with 
filters were used in what were essentially two separate but 
simultaneous scattering experiments on the same scat- 
tering volume. The incident and scattered beams were 
arranged such that the scattering angles 9b and 9g were 
different but the scattering vectors were identical, i.e. 
qs = Qg = q, where |q| = 4T7,7Tsin(6'/2)/A and n is the 
refractive index of cis-decalin. The outputs of the two 
detectors were cross-correlated to give the intensity cor- 
relation function (ICF): 



(JB(g,0)/G(g,r)) 
{lB{q)) {lG{q)) 



(1) 



where Ib and Iq are the scattered blue and green in- 
tensities respectively, and the angled brackets denote en- 
semble averages. 

In this arrangement, it can be shown [31] that only 
the singly-scattered light of each color probes exactly the 
same Fourier component of the density fluctuations and 
thus are correlated. All other, multiply-scattered, light 
does not probe the same component for both colors and is 
completely uncorrelated and thus does not contribute to 
the time-dependence of g^^\q,T). This can be expressed 
[31] in terms of the normalized single scattering dynamic 
structure factor, f(q,T): 

9^'\q,r) = l + (}'(}ls[fiq,r)f , (2) 

where the factor depends on the ratio of detector area 
and coherence area for single scattering and also on the 
overlap of the scattering volumes probed by each color. 
This factor is instrument related and dependent on scat- 
tering angle but not on the sample used. The other factor 
Pms reflects the fraction of singly-scattered intensities, 
{Ib) and (/q), relative to the total (singly and multi- 
ply) scattered intensities: 
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(4) (4) 

{Ib){Ig) 



(3) 



The concentration of polymer in our samples is low. The 
highest ratio of intensity scattered from polymer to that 
from colloid was measured to be 4 x 10~^. This high- 
est ratio only applied for one sample (H in Fig. 1) at 
the lowest scattering angle. Therefore we assume that 
the scattered intensity is from colloids only. Under these 
conditions, our measurements probe the structure and 
dynamics of the particles alone. 
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The static structure factor was measured with the pro- 
cedure described in [32]. First, the total average in- 
tensities, {Ib) and (Ig), and the intercept, 0^(3\is — 
(7^^^ ((7,0) — 1, of a concentrated sample of interest were 
measured at different scattering vectors q. The sample 
was rotated continuously during the measurement to en- 
sure ensemble average. Since the rotation only changes 
the time-dependence oi g'^^^q, r), the intercept and aver- 
age intensity wore not affected. It was found that it was 
necessary to average measurements at different height in 
the sample for non-ergodic samples to reduce random 
noise from the finite number of speckles sampled. Then 
the same measurements were made on a dilute suspen- 
sion of known (j) to obtain the single-particle form factor. 
The volume fraction of this dilute sample was (j)dn = 0.01, 
small enough that multiple scattering can be ignored, so 
that the measured intercept contains only the instrument 
related factor: = g^^^ (q, 0) — 1, which was the same as 
that in the measurement of the concentrated sample. 

The static structure factor is the ratio of singly- 
scattered intensity per particle from the concentrated 
sample to that from the dilute sample: S{q) = 
{{l^) /(j))/{{Idti) /4>dii)- This was calculated by taking 
into account multiple scattering and attenuation of light 
through the sample: 



S{q) = 
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where T is the transmission coefficient of the sample, and 

subscripts B, G arc for blue and green light respectively. 

Dynamic light scattering aims to measure the dynamic 
structure factor (DSF), /{q.r). This can be extracted 
from normalizing the ICF using Eq. 2: 



fiQ,r) 



/g(2)(g,T) - 1 



(5) 



However, since most of our samples were either non- 
ergodic or had very slow relaxation times, the time- 
averaged ICF only mcasurc;d fluctuations in a small sub- 
space of the whole configuration space. Explicit ensemble 
averaging was therefore required, and was performed by 
two methods. For short times (10^^ s < r < 20 s). 
brute- force ensemble averaging was done. Several hun- 
dred (typically between 500-865) of time-averaged ICF's, 

(2) 

gl {q,T), and associated scattered intensities, Ist and 
iGt, were measured, each for a duration of 40-60 s. Be- 
tween each measurement, the sample was rotated by a 
small angle to a different position so that each time- 
averaged ICF sampled a different Fourier component. 
The ensemble-averaged ICF was then constructed as: 



iBtiGtgt'^\q,T)^ 

(iBt) {iGt) • 

For longer times (r > 1 s), echo DLS was used. Details 

will be given elsewhere [33]. It essentially involves en- 
semble averaging by rotating the sample continuously at 



a constant speed and correlating the intensities at a small 
range of delay times around exact multiples, r ~ nT, of 
the rotation period, T, where the correlation function 
shows 'echo' peaks. The rotation decorrelates the ICF 
very quickly at small r. However, after a whole number 
of revolutions, the sample comes back to the same orien- 
tation and the ICF recovers to a value that is dependent 
only on the dynamics of the particles in the sample over 
that period of time. This gives rise to peaks in the ICF at 
r = nT. The maxima of these peaks follow the ensemble- 
averaged dynamics of the sample because the obtained 
ICF is an average over thousands of independent speck- 
les per revolution. We also correct for imperfections in 
the rotation using the area under each echo instead of 
the echo maximum [33]. The corrected ICF at r„ = nT 

was calculated from the measured ICF gin {q, t) as: 



S^'feO)-! 



+ 1, 



(7) 



where ^(t„) is the area under the echo around r„ = nT, 

A{t„) = J {gl^^ (g, r) - l)dr. The DSF was then obtained 
from corrected ICF in the usual way from Eq. 5. 

We used echo DLS to measure dynamics in the range 
r = 1 — 10^ seconds. Since the rotation used in- 
troduces slightly different alignment in the DLS setup 
(hence different /3^), the resulting intercepts are different 
from those obtained by brute-force ensemble averaging. 
Therefore we scaled the intercept of the echo DLS results 
by an arbitrary factor (in the range of 1-2) so that the 
resulting DSF from both methods matched in the region 
of overlap. 



III. RESULTS AND DISCUSSION 
A. Phcise diagram 

The equilibrium phase diagrams of colloid-polymer 
mixtures at different values of ^ are well known [9] . The 
non-equilibrium behavior of systems with ^ « 0.1 at 
low volume fractions {(p < 0.2) has been studied before 
[10, 11]. Here we concentrate on the higher volume frac- 
tion region (</> > 0.3). 

Many samples were prepared in a range of composi- 
tions of interest. After being homogenized by prolonged 
tumbling, samples were left undisturbed for observation. 
Because the sizes of colloidal particles are similar to wave- 
lengths of visible light, colloidal crystals can be seen with 
the naked eye as iridescent specks. 

Our observations are shown in Fig. 1. (These obser- 
vations have been presented and briefly discussed before 
[24].) In agreement with equilibrium theory [8] for sys- 
tems with short-range attraction, we observed an expan- 
sion of the fluid-crystal coexistence region upon increas- 
ing polymer concentration (diamonds). To the left of 
this region is a stable fluid phase (triangles) and to the 
right is the fully crystalline phase (inverted triangles). 
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These observations also agree with previous experiments 
on similar systems [9]. 

However, samples with very high colloid volume 
fractions and/or polymer concentrations (filled circles, 
squares and crosses) failed to crystallize for weeks to 
months even though equilibrium statistical mechanics 
predicts either fluid-crystal coexistence or full crys- 
tallinity. Samples with high colloid volume fractions and 
low polymer concentrations (circles in Fig. 1) showed 
all the characteristics of hard-sphere colloidal glass [2]. 
Weeks after homogenization and left undisturbed, sedi- 
mentation showed its effect: very thin layers at the top of 
the samples developed heterogeneous crystals due to the 
boundary effect of the meniscus and gravity. Samples de- 
noted by squares in Fig. 1 , with high polymer concentra- 
tion and moderate colloid volume fraction, showed signs 
of transient gels. They collapse under gravity after some 
'latency time' as observed previously in similar systems 
[10, 11]. However, the amount of collapse decreased and 
transient time increased dramatically in higher volume 
fraction samples. For concentrated samples with colloid 
volume fraction above 0.55, it took more than 4 weeks 
to see tiny collapses of less than half a millimeter at the 
very top of the meniscus. These collapses were distin- 
guished from normal sedimentation by their characteris- 
tic sharp and non-flat boundary between the collapsed 
material and a clear supernatant. No crystallization was 
observed in these samples however long they were left 
undisturbed. Interestingly, for non-crystallizing samples 
with very high colloid volume fraction and polymer con- 
centration (crosses), characteristics of both hard-sphere 
glass and transient gels were present. After 4-8 weeks, 
tiny collapses were seen, and also a thin layer of crystal 
phase appeared just under the collapsing boundary. 

Consider a sequence of samples of similar colloid vol- 
ume fraction and increasing attraction, for example sam- 
ples A-H in Fig. 1 with cf) ^ 0.6. According to thermody- 
namic equilibrium theory, all these samples should crys- 
tallize [8] . Sample A without polymer was a glass. Sam- 
ple B with a small amount of polymer was also a glass as 
no homogeneous crystallization was observed for 4 weeks 
and only hcterogeneously nucleated crystals at the menis- 
cus were observed after 13 days. However, sample C with 
~ 1.4 mg cm~^ of polymer completely crystallized in 1 
day. This means the glass transition line has moved to 
higher Failure to crystallize was seen again for sam- 
ples with polymer concentration above ~ 2.5 mg cm~^ 
(samples F,G,H). The behavior of all the samples in this 
region taken together show that the line of failure to crys- 
tallize had a re-entrant shape. 

In pure hard-spheres, crystallization ceases at essen- 
tially the same volume fraction as where /(g,oo) first 
becomes non-zero, i.e. at the glass transition [2, 3]. If 
this coincidence still holds for attractive hard-spheres 
systems, then we have observed a re-entrant glass tran- 
sition in hard spheres with short-range attraction. 

Previous studies of sticky hard spheres by MCT [16, 20] 
and computer simulation [17, 18] suggest that the re- 



entrant behavior is due to two different mechanisms of 
glassy arrest. The heuristic picture is as follows. In 
the 'repulsion-dominated' hard-sphere glass, particles are 
caged by their neighbors at high enough volume frac- 
tion. Short-range attraction clusters the particles of the 
cage and opens up holes, ultimately melting the glass. 
However, increasing the attraction further leads to an 
'attraction-dominated' glass where particles stick to their 
neighbors with long-lived bonds. In this terminology, 
samples A and B are repulsive glasses and F-H are attrac- 
tive glasses. Samples I-K must lie in the region where 
these two types of glass merge as they show characteris- 
tics of both types, with further evidence in the dynamics 
shown in section IIIC. The next sections, with results 
from light scattering, will give insights into the structure 
and dynamics of these glasses, and the nature of the re- 
entrant transition between them. 



B. Static structure factor 

We measured the static SSFs of the samples whose 
symbols are circled in Fig. 1. Note that samples C-E 
were measured as metastable fluids, i.e. before any crys- 
tal nucleation took place. Consider first the results for a 
sequence of samples (A-H) with (j) « 0.6, Fig. 2. These 
samples span the re-entrant glass transition line where 
the crystallization behavior showed dramatic changes. 
However, no re-entrant behavior can be seen in the SSF. 
Instead, there are only gradual changes upon increasing 
the attractive interaction. These gradual changes have 
been predicted by theory [16], and observed before in 
other experimental systems [34]. 

The most obvious and most easily quantifiable changes 
are in the height and position of the main peak. Broadly 
speaking, and taking experimental uncertainties into ac- 
count, the peak reduces in height and shifts to higher 
q when the attraction is increased (inset Fig. 2(a)). 
In detail, the peak position, q* , remains constant (at 
q* R ~ 3.8; samples A D) until just before we enter the 
attractive glass region (sample E), whereupon it increases 
by « 5% to reach another constant value {q*R ~ 4, sam- 
ples F-H). These samples have approximately constant 
(j) (in fact it decreases slightly from A to H, Fig. 1). The 
increase in q* is the result of a significant fraction of 
neighboring particles becoming trapped in each others' 
narrow depletion potential well when the attractive glass 
forms. Quantitatively, a 5% increase in q* corresponds to 
a 15% increase in the local packing fraction, from 0.6 to 
0.69; the latter is the random close packing volume frac- 
tion for our system (measured by spinning down a sample 
of known </>). In other words, the nearest particles in the 
attractive glass are practically touching. 

The clustering of particles at constant volume necessar- 
ily implies that the average number of nearest neighbors 
should decrease, and that 'holes' are opened up to render 
the structure more inhomogeneous on the spatial scale of 
a few particles. The former is reflected in the decrease in 



5 



S{q*). Significantly, upon increasing the attraction from 
zero, the decrease in the peak height starts at the point 
of the melting of the repulsive glass, and continues until 
we enter the attractive glass region, whereupon the peak 
height remains constant (inset Fig. 2(a)). 

The increased heterogeneity is reflected in a rise in the 
SSF at low q, Fig. 2(b). The smallest q we have stud- 
ied was qminR = 1.50, corresponding to a length scale 
of about 4 particle radii. The value of S{qmin) increases 
nearly exponentially with the polymer concentration be- 
tween samples A-E (inset Fig. 2(b)), and thereafter re- 
mains constant. The increased density fluctuations at 
this length scale corresponds to the opening up of 'holes' 
due to particle clustering. 

Note that all three features considered, q* , S{q*) and 
S{qmin), remain virtually constant for all three attractive 
glass samples, F-H. Once particles drop into each others' 
narrow attractive potential wells, any further structural 
changes will be hard to resolve. We shall see, however, 
that the dynamics continues measurably to evolve from 
sample F to sample H: in this regime of almost-touching 
nearest neighbors, a very small change in the structure 
has very large dynamic consequences. 

All the qualitative features we observed in the evolu- 
tion of the SSFs for samples A-H are also seen in the 
SSFs for samples I-K at the higher volume fraction of 
(j) « 0.64, Fig. 3. However, the effects are significantly 
less obvious, largely because the range of polymer con- 
centration is now much smaller and (j) is higher. At low q, 
the values of /S'(gmin) are lower than those of similar poly- 
mer concentration but lower (j> (C-E) (Fig. 2(b)). This is 
because at higher volume fraction, a tight local clustering 
of some particles does not create so much room elsewhere 
— there is less space for developing heterogeneities. 



C. Dynamic structure factor 

Our goal is to study how the polymer-induced deple- 
tion attraction affects the particle dynamics. But the 
presence of the polymer influences the dynamics in an- 
other, essentially trivial, manner — by increasing the ef- 
fective viscosity of the medium in which the particles dif- 
fuse from that of the pure solvent, rjo, to that of a polymer 
solution, r]rr]o at concentration Cp**®. To determine r]r, we 
measured the viscosity of pure polymer solutions with a 
miniature suspended-level viscometer. Fig. 4, and used 
a quadratic fit to the data to obtain rjr for our samples. 
The value of ''''in each sample was estimated from Cp, (j) 
and ^ using an approximate expression based on scaled- 
particle theory [8]. 

The rate of dynamical decay at wave vector q depends 
on the length scale being probed; in dilute systems it 
scales as q^ . Thus, in order to compare the dynamics of 
different samples at different wave vectors, and to high- 
light the effects of the attraction, we scaled the delay 
time variable of the DSFs by the relative viscosity rjr and 
the dimensionless wave vector {qR}^, so that DSF is pre- 



sented as a function of the 'scaled time' {qRYr/rir. Note 
that for the lowest q studied, the scaled time is very close 
to the real time, while at the highest q, it is increased by 
about an order of magnitude. 

We found aging [35] in all non-crystallizing samples. 
The dynamics slowed down with the 'waiting time' — 
the time interval between the cessation of tumbling and 
the beginning of measurements. Fig. 5. It is known that 
the hard-sphere glass ages [36] . We found that the rate of 
aging in different glasses were different and that its effects 
were complex. Repulsive glasses aged only in the first day 
or two, after which they did not evolve within the time 
window of the measurements. Attractive glasses, on the 
other hand, showed different dynamics with age for up to 
10 days. Aging is complicated enoTigh to be the subject 
of a separate study and was not investigated systemat- 
ically in this work. To eliminate as much as possible 
aging effects on dynamical results within practical limits 
of waiting time, we present DSFs of glassy samples with 
age betwecm 1 and 4 days. The dynamics of crystalliz- 
ing samples (C-E) were measured while they were in the 
metastable state well before the appearance of crystal- 
lization. Below we first show results of different samples 
at the same q, then at different q for the same sample. 



1. Constant scattering vectors, variable compositions 

The DSFs of samples A H at gi? = 1.50, Fig. 6, clearly 
evolve non-monotonically with increasing polymer con- 
centration and show re-entrant behavior. Briefly, sam- 
ples A and B are non-ergodic within our time window, 
while samples C-E are ergodic (their DSFs decay com- 
pletely to zero), and samples F-G become non-ergodic 
again. 

In detail, the DSF of sample A, a pure hard-sphere 
glass, shows a plateau at /^(g, oo) w 0.7, correspond- 
ing to particles getting 'stuck' in their nearest-neighbor 
cages. This can be compared with previous work [3, 37]. 
Note that in doing so, it is important to compare samples 
with the same density relative to random close packing: 
i.e. the same {(j)rcp — (f>)/<Picp, since (prep differs according 
to the polydispersity of the colloids [38]. 

With a small amount of polymer added to the hard- 
sphere glass, sample B shows the same qualitative dy- 
namics. Quantitatively, however, the height of the 
plateau is lower, /(g, oo) — 0.62. This indicates that 
particles in B are not as restricted as in A, i.e. the cage 
is loosened by the attractive interaction, but still remains 
closed in our time window. 

The DSF of sample C decayed completely in (a 'scaled 
time' of) about 1000 seconds, as did those for the other 
crystallizing samples D and E. It is interesting to note 
that the DSFs of these three samples slow down upon 
increasing polymer concentration but all reach zero at 
about the same scaled time. The DSF of sample C shows 
the remnant of a plateau at a scaled time of « 10 s. The 
DSFs for samples D and E exhibit a very stretched single 
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decay, rather than a two-stepped process. This is unusual 
behavior for a fluid at volume fraction </> ~ 0.6 (at least 
at first sight). 

The intermediate, /?, and long-time, a, decay in a dense 
hard-sphere fluid are attributed to particles 'rattling' in 
their local neighbor cages, and escaping from these cages, 
respectively [3] . Attraction hinders the 'rattling' by trap- 
ping particles in potential wells, but accelerates the cage 
opening by clustering. At some polymer concentration 
(or attraction strength), the two time scales coincide. If 
at this point the attraction alone is not enough to trap 
the system in an non-ergodic state, wc will observe the 
melting of the repulsive glass into an ergodic fluid domi- 
nated by attraction. This is the case for sample C, where 
the a and /3 decays are barely distinguishable in the DSF. 
At higher polymer concentrations, the cage concept is no 
longer appropriate for describing the particle dynamics 
— for it to be valid, a particle has to 'rattle' many times 
in a cage before it opens. Note that this is a distinctive 
feature of dense fluids with short-range attraction. In a 
dense fluid with long-range attraction, the effective po- 
tential well experienced by any particle due to its neigh- 
bors is essentially flat. This adds a (negative) constant 
to the free energy, so the phase behavior [39] and dy- 
namics of the system are still controlled by repulsion (or, 
equivalently, entropy) . 

Note that the shape of the DSF of sample C at gi? = 
1.50 is similar to that shown in curve 2, Fig. 11 of [16]. 
This DSF was calculated at qR = 2.1 for a sample in 
the re-entrant portion of the state diagram in a system 
that just shows a glass-glass transition and an A3 point. 
Recent calculations for colloid-polymer mixtures [20] sug- 
gests that our system, with ^ w 0.09, should show exactly 
these features. 

The DSFs of samples F-H are, once more, non-ergodic 
in our time window: they do not decay completely even 
after 10** seconds. Simple extrapolation indicates that it 
would take these DSFs at least 10^ seconds to reach zero. 
The DSFs of samples G and H show points of inflection; 
that for sample H is clearer and occurs at / = 0.995 — a 
very high value compared to the plateaus in hard-sphere 
glasses. These high points of inflection can be associated 
with dynamics originating from particles rattling in very 
narrow attractive potential wells. 

At other wave vectors. Figs. 7 and 8, the DSFs behave 
in a similar way, namely relatively low plateaus in the re- 
pulsive glasses A and B, complete decay in the metastable 
fluids C-E, and extremely slow dynamics in the attrac- 
tive glasses F-H, with very high points of inflection in 
G and H. Note, however, that at the peak of the cor- 
responding SSFs, the DSFs for samples C-E are barely 
distinguishable (Fig. 7). 

The plateaus in the DSFs of the repulsive glasses can 
be used as a measure of f{q, oo), the non-ergodicity pa- 
rameter. An estimate of this quantity for the attractive 
glasses is more problematic, partly due to significant ag- 
ing in our time window. To proceed, wc use the value of 
/ at the point of inflection as a surrogate; we call this 



the 'measured' f{q, oo) = f^^'^^q, oo). The evolution of 
y(*f)(^g^oo) with increasing polymer concentration (sam- 
ples A-H) is shown in Fig. 9. The non-ergodicity pa- 
rameter decreases slightly when moving from A to B, 
away from the hard-sphere glass. When attraction melts 
the repulsive glass, f^^\q, oo) = for samples C-E 
(not shown). Sample F did not crystallize and showed 
non-ergodic dynamics up to 10* seconds but did not ex- 
hibit any discernible point of inflection in its dynamics. 
Samples G and H had extremely high non-ergodicity pa- 
rameters of nearly 1. A 'jump' in f{q, oo) when moving 
from repulsive to attractive glass was predicted by MCT 
(Fig. 7 in [16]). 

The evolution of the short-time dynamics of the whole 
sequence of samples is also interesting. Fig. 10 shows 
the short-time behavior of the DSFs for A-H at large 
length scale, qR = 1.50, where experimental noise is low- 
est. Note the very small vertical interval, 1.000 to 0.997, 
spanned in this figure; thus only the first 0.3% of the 
decays of the DSFs are being analyzed. The DSFs of re- 
pulsive glasses A and B possessed relatively long linear 
parts, corresponding to the first term in r in the ex- 
pression derived from the Smoluchowski (many-particle 
diffusion) equation [40]: f{q,T) = ^ji^q^r + Oir^), 
where Dq is the free-particle diffusion constant in pure 
solvent (with no polymer) Dq — ksT/GwrjoR, and H{q) 
is the hydrodynamic factor. This linear regime of the 
DSFs indicates that at short time, individual particles 
still diffuse freely without the influence from direct in- 
teraction with their neighbors. The change in limiting 
slope as T ^ 0, or the short-time diffusion coefiicient 
Ds{q) = DoH{q)/ S{q), can be almost entirely explained 
by the change in S{q) (Fig. 11), including the strong de- 
crease on entering the attractive glass regime. What is 
more interesting is that the dynamics depart from free 
diffusion progressively earlier upon increasing attraction 
(Fig. 10). In fact, for the attractive glasses F H. the par- 
ticles are confined so tightly by the attractive potential 
wells that the DSFs display non-linearity almost imme- 
diately (cf. also insets to Figs. 6-8). 

Moving to the (shorter) sequence of samples at the 
higher volume fraction of (j) w 0.64 and closer to the in- 
tersection of the two glass transition lines, samples I -K 
in Fig. 1, we see the emergence of remarkably stretched- 
out, extremely slow dynamics. Consider first the data at 
qR = 1.50, Fig. 12. In terms of short-time dynamics (in- 
set, Fig. 12), samples I and J arc comparable to samples 
C and D, while sample K shows a behavior intermediate 
between those of samples E and F. At intermediate times, 
the decay is linear with respect to the logarithm of the 
scaled time. Thereafter there is an incipient plateau at 
/ 0.7 in sample I, reminiscent of the plateau in repul- 
sive glasses A and B, before a further decay, but never 
beyond ^ 0.62 in our time window. There is no incipient 
plateau for the other two samples. Note that the DSF 
of sample I shows aspects of the behavior of repulsion- 
dominated glasses (long time) and a fluid dominated by 
short-range attraction (short time) . The two regimes are 
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'bridged' by a stretched log-time decay. 

At the peak of the SSF, Fig. 13, sample I behaves in 
a similar way at short to intermediate times, while there 
is no incipient plateau at long times. Samples J and K 
develop an incipient plateau as high as ~ 0.993 (inset 
Fig. 13) before turning over to decay more rapidly in 
logarithmic time. 

The fact that these samples show extremely stretched 
out dynamics, logarithmic in time, suggests that they 
are very close to the A3 critical point predicted by MCT, 
where the repulsive and attractive glasses become indis- 
tinguishable [14, 16, 21]. In particular, the shape of the 
DSF of sample I at gi? = 1.5 is comparable to curve 3 
in Fig. 11 of [16], calculated at qR = 2.1 for a sample on 
the repulsive glass transition line very close to where it 
intersects the attractive glass transition line for a system 
that just shows an A3 singularity. This is not inconsis- 
tent with the position of sample I on the state diagram. 
Fig. 1, of our system at ^ « 0.09 [20]. 

Heuristically, we may begin to make sense of log-time 
decays as follows. At high enough volume fraction, the 
average distance between neighboring particles will de- 
crease to a value such that they are always well within 
the attraction range of each other [41]. If the attraction 
is strong enough, the restriction of particle movement 
due to the neighbor cage and the restriction caused by 
bonding between particles take place simultaneously at 
all times. This competition between two opposite mech- 
anisms may lead to a broad distribution of decay times 
and therefore a very stretched out DSF [42] . 



2. Constant compositions, variable scattering vectors 

In this section, we show for completeness the dynamics 
of each sample at different scattering vectors in Figs. 14- 
17. The change of DSFs with q in repulsive glasses A and 
B are in agreement with previous work [3, 37]. Other 
samples show the general trend that the dynamics be- 
come slower at scattering vectors with higher S{q). The 
only exception concerns the intermediate-time dynamics 
of the attractive glasses F-H (insets, Fig. 16). The sig- 
nificance of the rather complicated (jr-dependcncc of the 
intermediate-time dynamics of these samples is not clear. 
Nor do we know of any detailed calculations to date that 
can throw light on this issue. 

The systematic g-dependent data shown in Figs. 14- 
17 allow us to investigate the ^'-dependence of the mea- 
sured non-ergodicity parameter, f^^\q, oo), in detail. 
The measured non-ergodicity parameters of glassy sam- 
ples A, B, G and H are shown as a function of scattering 
vector q in Fig. 18. The data for repulsive glasses A 
and B vary essentially with the static structure factor, 
as observed in hard-sphere glasses [37]. Attractive glass 
G and H on the other hand showed extremely high mea- 
sured non-ergodicity parameters that hardly vary with q. 
This agrees with predictions by MCT (c.f. Fig. 8 in [16]). 



IV. CONCLUSION 

We have studied a dense system of hard-sphere col- 
loids with a short-range inter-particle attraction induced 
by the depletion effect of added non-adsorbing polymer. 
The observed crystallization behavior as well as parti- 
cle dynamics stiidied by DLS reveal a re-entrant glass 
transition. With little attraction, the system at high 
enough volume fraction is 'stuck' in a repulsive glassy 
state where the arrest is due to caging by neighboring 
particles. Our data support the suggestion [13] that at- 
traction causes particles to cluster, thus opening up holes 
in the cages and melting the glass. At the same time, 
the attraction slows down the particle dynamics. We 
found that the repulsive glass melts when the character- 
istic time of the attraction-dominated particle dynamics 
becomes comparable to that of cage opening. The result- 
ing ergodic fluid shows a distinctive dynamical feature: 
despite the fluid's high density, its DSF does not show 
distinct a and P relaxation processes. Increasing the at- 
traction further leads to different kind of arrest where 
the strong attraction between particles create long-lived 
bonds and prevent structural rearrangement, giving rise 
to an attraction-dominated glass. Detailed light scatter- 
ing has been used to probe the effect of attraction on 
both structure and dynamics. 

Qualitatively, this scenario agrees well with predictions 
from MCT calculations (with those reported in [20] be- 
ing closest to the present experimental system). In par- 
ticular, we observed very slow, log-time dynamics in the 
DSFs in the region where the two glass transition lines 
are expected to meet. Quantitatively, however, our re- 
sults stand as a challenge to MCT (or any other theory): 
the detailed calculations needed for direct quantitative 
comparison have not, to our knowledge, been performed. 

A detailed comparison between experiment and the- 
ory faces a number of non-trivial problems. First and 
foremost, since calculated and measured glass transition 
thresholds differ, choices exist as to what constitute 'cor- 
responding state points' for the purpose of making the 
comparison. In the case of pure hard spheres, where 
^MCT _ 0.52 and (f)^''^^ « 0.58, it is accepted practice 
to compare measurements and calculations at the same 
relative volume fraction {(p — (pg)/<pg [3]. The situation is 
more complex in a colloid-polymer mixture, since a state 
point is now specified by the densities of both compo- 
nents. The predicted glass transition lines show quan- 
titative disagreement with experiments over the whole 
composition plane (cf. Fig. 1 in [17]). To compare cal- 
culated and measured SSFs and DSFs, a protocol for 
identifying 'corresponding state points' is needed. 

Secondly, the attractive interaction between two par- 
ticles is always specified directly as a potential energy in 
calculations. The corresponding experimental variable 
is the polymer concentration in the free volume, c^'^''. 
This is currently guessed at using an uncontrolled and 
untested approximation based on scaled-particle theory 
[8] , and is likely to lead to large systematic errors in sys- 



8 



tems with high coUoid volume fractions. Thirdly, the 
marked and complex aging behavior of the attractive 
glasses complicates the definition of a non-ergodic state 
for the purposes of comparing with MCT. Despite these 
potential difficulties, however, our data suggest that it 
may be worthwhile performing a series of calculations at 
fixed (j) and increasing attraction crossing the re-entrant 
gap in between the repulsive and attractive glass transi- 
tion lines for a system of hard spheres interacting with 
something like an Asakura-Oosawa potential [20]. 

Finally, it is clear that attractive and repulsive glasses 
show qualitative distinct aging behavior. Classical MCT 
does not predict aging, but it is a generic feature of ex- 
perimental glasses of all kinds [35] . A number of theoret- 
ical approaches are emerging, and simulation is a valu- 
able tool. It is probable that further study of this phe- 



nomenon in our model colloid-polymer mixture should 
throw significant light on this intriguing (and generic) 
phenomenon [22]. 



Acknowledgments 

KNP is funded by an UK ORS award and the Univer- 
sity of Edinburgh. Partial financial support for experi- 
mental equipment came from EPSRC grant GR/M92560. 
We thank Matthias Fuchs, Michael Cates, Antonio Puer- 
tas, Johan Bergenholtz and Francesco Sciortino for il- 
luminating discussions during various stages of this 
work, and Michael Cates also for commenting on the 
manuscript. 



[1] C. A. Angell, Science 267, 1924 (1995). 
[2] P. N. Pusoy and W. van Mogcn, Nature 320, 340 (1986). 
[3] W. van Megon and S. M. Underwood, Phys. Rev. E 49, 
4206 (1994). 

[4] U. Bengtzelius, W. Gotze, and A. Sjolander, J. Phys. C 

17, 5915 (1984). 
[5] W. Gotze and L. Sjogren, Rep. Prog. Phys. 55, 241 

(1992). 

[6] D. Rudhardt, C. Bechinger, and P. Leiderer, Phys. Rev. 

Lett. 81, 1330 (1998). 
[7] S. Asakura and F. Oosawa, J. Chem. Phys. 22, 1255 

(1954). 

[8] H. N. W. Lekkerkerker, W. C. K. Poon, P. N. Pusey, 
A. Stroobants, and P. B. Warren, Europhys. Lett. 20, 

559 (1992). 

[9] S. M. Ilett, A. Orrock, W. C. K. Poon, and P. N. Pusey, 

Phys. Rev. E 51, 1344 (1995). 
[10] W. C. K. Poon, A. D. Pirie, and P. N. Pusey, Faraday 

Discuss. 101, 65 (1995). 
[11] W. C. K. Poon, L. Starrs, S. P. Meeker, A. Moussaid, 

R. M. L. Evans, P. N. Pusey, and M. M. Robins, Faraday 

Discuss. 112, 143 (1999). 
[12] J. Bergenholtz and M. Fuchs, J. Phys.: Condens. Matter 

11, 10171 (1999). 
[13] J. Bergenholtz and M. Fuchs, Phys. Rev. E 59, 5706 

(1999) . 

[14] L. Fabbian, W. Gotze, F. Sciortino, P. Tartaglia, and 
F. Thiery, Phys. Rev. E 59, R1347 (1999). 

[15] F. Mallamace, P. Gambadauro, N. Micali, P. Tartaglia, 
C. Liao, and S.-H. Chen, Phys. Rev. Lett. 84, 5431 

(2000) . 

[16] K. Dawson, G. Foffi, M. Fuchs, W. Gotze, F. Sciortino, 

M. Spcrl, P. Tartaglia, T. Voigtmann, and E. Zaccarelli, 

Phys. Rev. E 63, 011401 (2001). 
[17] K. N. Pham, A. M. Puertas, J. Bergenhohz, S. U. Egcl- 

haaf, A. Moussaid, P. N. Pusey, A. B. Schofield, M. E. 

Cates, M. Fuchs, and W. C. K. Poon, Science 296, 104 

(2002). 

[18] A. Puertas, M. Fuchs, and M. E. Gates, Phys. Rev. Lett. 

88, 098301 (2002). 
[19] T. Eckert and E. Bartsch, Phys. Rev. Lett. 89, 125701 

(2002). 



[20] J. Bergenholtz, W. C. K. Poon, and M. Fuchs, Langmuir 

19, 4493 (2003). 
[21] F. Sciortino, P. Targaglia, and E. Zacarelli, preprint, 

cond-mat/0304100. 
[22] E. ZacarelU, G. Foffi, F. Sciortino, and P. Tartagha, 

preprint, cond- mat/0304192. 
[23] W. C. K. Poon, J. Phys.: Condens. Matter 14, R859 

(2002) . 

[24] W. C. K. Poon, K. N. Pham, S. U. Egelhaaf, and P. N. 

Pusey, J. Phys.: Condens. Matter 15, S269 (2003). 
[25] L. Antl, J. W. Goodwin, R. D. Hill, R. H. Ottewill, S. M. 

Owens, S. Papworth, and J. W. Waters, Colloids Surf. 

17, 67 (1986). 

[26] S. Auer, W. C. K. Poon, and D. Frenkel, Phys. Rev. E 

67, 020401 (2003). 
[27] P. N. Pusey and W. van Megen, J. Chem. Phys. 80, 3513 

(1984). 

[28] G. C. Berry, J. Chem. Phys. 44, 4550 (1966). 

[29] P. G. Bolhuis and D. A. Kofke, Phys. Rev. E 54, 634 

(1996). 

[30] P. Bartlott, J. Chem. Phys. 107, 188 (1997). 

[31] P. N. Segre, W. van Megen, P. N. Pusey, K. Schatzel, 

and W. Peters, J. Mod. Optics 42, 1929 (1995). 
[32] A. Moussaid and P. N. Pusey, Phys. Rev. E 60, 5670 

(1999). 

[33] K. N. Pham, S. U. Egelhaaf, A. Moussaid, and P. N. 

Pusey, to be published. 
[34] T. Eckert and E. Bartsch, Faraday Discuss. 123, 51 

(2003) . 

[35] J.-P. Bouchaud, in Soft and Fragile Matter and Nonequi- 

librium Dynamics and Metastability and Flows, edited 
by M. E. Cates and M. R. Evans (Institute of Physics 
Publishing, 2000), Scottish Universities Summer School 
in Physics. 

[36] W. van Megen, T. C. Mortensen, S. R. WiUiams, and 

J. MiiUer, Phys. Rev. E 58, 6073 (1998). 
[37] W. van Megen, S. M. Underwood, and P. N. Pusoy, Phys. 

Rev. Lett. 67, 1586 (1991). 
[38] W. Schaort and H. Sillescu, J. Stat. Phys. 77, 1007 

(1994). 

[39] B. Widom, Science 157, 375 (1967). 

[40] P. N. Pusey, in Liquids, Freezing and the Glass Transi- 



tion, edited by J. -P. Hansen, D. Levesque, and J. Zinn- 
Justin (Elsevier, Amsterdam, 1991), chap. 10. 

[41] In our system the estimated distance between particles 
for samples I-K (0 — 0.64) from random close pack 
{(prep = 0.69) is (0rcp/</>)^'^^ = 1-03 ~ 1 + ^/S, where 
the attractive potential is half of the maximum depth. 

[42] Formally, a distribution of decay times gives a decay 
linear in logr. Limitations in our data mean, however, 
that we cannot use them to back out the actual decay- 
time distribution in our samples. 



E 



A 


stable fluid 


O 


fluid-crystal 


V 


fully crystallised 


■ 


gel/attractive glass 


+ 


(both) glass 


• 


repulsive glass 



8)H 
FH 



/ ®B 



A oo v*M • mA 

^ \ , , , , \ , L. 



0.4 



0.5 

volume fraction ( 



0.6 



0.7 



FIG. 1: Equilibrium and non-equilibrium behaviors of a 
colloid-polymer mixture of ^ = 0.09. Open symbols are those 
that reached thermal equilibrium (fluid, fluid-crystal coexis- 
tence, and fully crystallized). Other samples did not crystal- 
lize: some showed characteristics of hard-sphere glasses at the 
onset of sedimentation (filled circles), some showed those of 
attraction-driven glasses and gels (filled squares), and some 
showed both (pluses). Dashed curves are guides to the eye 
showing the observed boundary where crystallization ceased. 
Light scattering data for marked samples labeled A-K are 
shown in the following figures. This diagram has been shown 
in [24]. 
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FIG. 2: Static structures factor of samples A-H {(p ~ 0.6) as 
function of dimensionless wave vector qR. The lines are guides 
to the eye. (a) The peak position q* shifts to higher q, while its 
height reduces and width increases upon increasing attraction. 
The inset shows the peak positions and heights of these static 
structure factors as a function of polymer concentration, (b) 
The same SSFs plotted with logarithmic vertical axis shows 
the increase of S{q) at low q. The inset shows S{q) at the 
lowest wave vector qR — 1.50. Vertical dashed lines in both 
insets indicate the glass transitions observed in Fig. 1. 
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FIG. 3: Static structure factors of samples I-K with (j> ~ 0.64. 
The inset shows the same data with a logarithmic vertical 
axis. 
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FIG. 4: The viscosity of polystyrene in cis-decalin at 20° C at 
different polymer concentrations. A quadratic fit (solid line) 
was used to interpolate to viscosities of samples with diflferent 
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FIG. 5: Aging in samples A, G and ¥L at qR ^ 2.93. The 
legend indicates the waiting time between the end of tumbling 
and the start of measurements. The DSFs slow down and the 
points of inflection become clearer with increasing age of the 
samples. 
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FIG. 6: Collective dynamic structure factors at qR — 1.50 
from samples A-H spanning the re-entrant region. The time 
axis is scaled to dimensionless length scale {qRf and relative 
polymer solution viscosity rjr. The inset shows the same plots 
on an expanded vertical axis. 
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FIG. 7: DSFs at the peak of the SSFs for samples A-H. The 
inset shows the same plots on an expanded vertical axis. 
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FIG. 8: DSFs of samples A-H at qR = 4.30, to the right of all 
S{q) peaks. The inset shows the same plots on an expanded 
vertical axis. 
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FIG. 9: The measured non-ergodicity parameters, 
/(^^(g, oo), at different wave vectors as a function of 
polymer concentration in samples (left to right) A, B, G, and 
H. The dashed lines indicate the glass transitions observed 
in Fig. 1. 




FIG. 10: The short-time dynamics of samples A-H at qR = 
1.50. The straight hnes are fits to the hnear part of the DSFs 
at r ^ 0. The dynamics departs from an initial diffusive 
regime progressively earlier upon increasing attraction. The 
short-time diffusion coefficient in the limit r ^ is also re- 
duced significantly. 
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FIG. 11: The normalized short-time diffusion coefficient 
Doze's (circles, left scale), static structure factor S[q) and hy- 
drodynamic factor H{q) (right scale) at qR = 1.50. Ds and 
S(g) were extracted from Fig. 10 and 2 respectively. The de- 
crease in Ds is nearly in line with the increase in S{q) so that 
H[q) only increased slightly. 
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FIG. 12: The DSFs at qR = 1.50 for samples I-K with 
(j) ~ 0.64. Extremely stretched relaxation is found in all three 
samples with logarithmic decay over long ranges of r (straight 
lines). The inset shows the short-time dynamics, which devi- 
ate from the diffusive regime from very early times. 
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FIG. 13: DSFs at the peak of the SSH for samples I-K. AU 
decay much slower than at low q. Sample I shows a logarith- 
mic decay for about 3 decades in scaled time. Samples J and 
K develop very high plateaus (inset). 
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FIG. 14: The DSFs of samples A and B at different q vectors, 
q* denotes the peak position of the static structure factor 
^(g). The general shape of the DSFs are very similar. Hori- 
zontal lines denotes the height of the plauteau (non-ergodicity 
parameters) that are plotted in Figs. 9 and 18. 
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FIG. 15; The DSFs of samples C-E at different q. The rate 
of decay varies in the opposite direction to S{q) (c.f. Fig. 2). 
However, all decay to zero at approximately the same scaled 
time. Except for sample C at the lowest q, all other DSFs 
do not show two distinct relaxation processes as other dense 
fluids. 
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FIG. 16: The DSFs of samples F-H at different q. The vertical 
axes span different ranges. Sample F did not show a point of 
inflection, but G and H have very high points of inflection 
(horizontal lines), the values of which are used in Figs. 9 and 
18. 
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FIG. 17: The DSFs of samples I, J and K at different q. The 
insets show the same quantities with expanded vertical axes. 
The relaxations show similar behavior at all wave vectors ex- 
cept at the peak of S{q). Sample I decays to a logarithmic 
section and then appears to turn up to a plateau. Sample J 
shows a very long section of logarithmic decay. Sample K is 
similar to J with a shorter stretch of logarithmic decay. In 
the early decay at the peak of S{q), the DSF of sample I has 
a long stretch of logarithmic decay whereas samples J and K 
develop very high plateaus. 
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FIG. 18: The measured non-ergodicity parameters of samples 
A, B, G and H as a function of scattering vector q (points), 
and the static structure factor of sample B (line) for compar- 
ison. The non-ergodicity parameters of repulsive glasses A 
and B follow the static structure factor, whereas those of the 
attractive glass are extremely high and hardly fluctuate with 
q (upper pane with expanded vertical axis). 



